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Abstract 

The value of a s (Mz) emerging from the so called global fits based mainly on 
the data at the Z peak (and assuming the standard model) is three standard de- 
viations higher than the one stemming from the low-energy phenomenology. The 
corresponding value of Aq CD is very large, ~ 500 MeV, and is incompatible with 
crucial features of QCD. If persists, the discrepancy should be interpreted as due 
to contributions to the Z-quark-antiquark vertices which go beyond the standard 
model. 



The statement that precision measurements of electroweak physics at the Z peak 
give no evidence whatsoever of new physics lying beyond the standard model (SM) is 
becoming common place now. Moreover, the values of the SM parameters extracted 
from analysis of the LEP and SLD data are used as canonical. This refers, in 
particular, to the strong coupling constant. The so called global fits which assume 
validity of the standard model and are based on a large set of data ("high-energy 
data") yield values of a s (M z ) in the MS scheme which cluster around 0.125 [p], U), 
with the error bars 0.005 ||. The corresponding value of Aqcd is about 500 MeV 0. 
These numbers, accepted as the most exact results for the strong coupling constant 
existing at present, propagate further into a stream of papers, published in the last 
year or two, devoted to various aspects of QCD. 

The question arises whether Quantum Chromodynamics can tolerate these num- 
bers. I will argue below that the answer is negative. There are two reasons why I 
believe that a s (M z ) must be close to 0.11 and the corresponding value of Aq CD close 
to 200 MeV (or even smaller). A more formal argument comes from consideration 
of traditional "low- energy" data, the cleanest of which is the evolution of moments 
of the structure functions in deep inelastic scattering (DIS). These measurements 
are abundant and have high statistics (for a review see ||). Theoretical formulae 
for the moment evolution are known in the leading and next-to-leading logarithmic 
approximations; if one considers the Euclidean domain of Q 2 above ~10 GeV 2 non- 
perturbative effects play no role in the Q 2 evolution of the moments provided that 
the moments considered are not too high. A typical result for a s emerging in this 
way (scaled to the normalization point /i = Mz) is 0.113 ± 0.005 0. 

Other "low-energy" analyses - heavy quarkonia, jets at PEP and PETRA and 
so on - produce similar and even lower values of a s , with the only exception to be 
discussed in some detail below. However, in determining the value of a s from these 
data one encounters serious problems: due to the essentially Minkowskean nature of 
the jet and quarkonium calculations precise estimates of the role of non-perturbative 
effects are difficult and, hence, the corresponding results for a s are inconclusive. 
Moreover, in some instances above even the perturbative next-to- leading corrections 
are not calculated so far. It should be noted that lattice calculations (e.g. Ref. ||) 
also produce a s {Mz) = 0.115 or less; being a lattice outsider I do not know how 
reliable these lattice results are. 

The second, less formal, argument in favor of Aq CD ^ 200 MeV is the success 
of the operator product expansion (OPE) || in an extremely wide range of ap- 
plications to different QCD problems. Although less formal at the moment, this 
argument seems more convincing to me. The analysis is based on the Euclidean 
expansions (QCD sum rules [[UJ) ensuring, thus, the best possible control over non- 
perturbative contributions. Although Wilsonian OPE is valid in any consistent field 
theory, practically successful numerical predictions become possible only because of 
the existence of a window in QCD - a crucial phenomenon not completely understood 
theoretically |fLl|| . The window is a Euclidean domain of momenta where in certain 



correlation functions perturbative corrections turn out to be numerically smaller 
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than non-perturbative ones. In typical instances after the Borel transformation the 
window stretches down to ~ 0.7 GeV. This fact ensures a very fast transition from es- 
sentially perturbative regime to essentially non-perturbative one. An interpretation 
of the window phenomenon emerges if one assumes that b, the first coefficient in the 
Gell-Mann-Low function, is a large parameter (numerically large), all perturbative 
corrections are suppressed by powers of 1/6, while the non-perturbative corrections 
do not contain this parameter. The window phenomenon would be impossible if 
Aq CD > 500 MeV. It should be stressed that the difference between a s (M z ) = 0.125 
and a s (Mz) = 0.11 (i.e. Aqcd ~ 500 MeV versus 200 MeV) is not merely quantita- 
tive but, rather, qualitative. In the first case one expects a rather slow and gradual 
transition from the perturbative regime to the non-perturbative one, while in the 
second case the non-perturbative (power) corrections blow up in the domain where 
perturbation theory still seems convergent. 

Precise measurements of Aqcd from the QCD sum rules can be fully formalized. 
As a matter of fact, the old work [|l!| in this direction is seminal. It yields Aqcd < 
210 MeV at the one-loop level. Later analyses along these lines also exist. Further 
efforts, both theoretical and experimental, are needed in order to measure Aqcd i n 
this way at the level of accuracy desirable today. But even the achieved level of 
accuracy rules out A QC d ^ 500 MeV. 

A discrepancy between the low-energy expectations for a s and the fits at Z 
alerted some theorists a few years ago (see below). Their arguments were largely 
overshadowed later by the assertion, worked out in a series of interesting and stimu- 
lating papers |L3]], that a precise low-energy determination of a s from r is possible, 
a s (M T ) = 0.33 ± 0.03. This result implies, in turn, that a s (M z ) = 0.120 ± 0.003, 
in accord with the value measured at the Z peak. Ref. [|l3j presents the state of 
the art in estimating all known sources of non-perturbative corrections. The leading 
corrections come from the gluon and four-quark condensates and turn out to be 
negligibly small. This calculation, quite correct by itself, gives rise to a new doc- 
trine - a very precise determination of a s (Mz) from essentially perturbative formula 
for r(r — > hadrons) is possible - and we are witnessing now how this doctrine is 
gradually becoming generally accepted in the community. 

The problem with all approaches of this type is that it can not be formulated as a 
completely Euclidean analysis. One has to deal with Minkowskean spectral densities 
integrated with some weights over some finite energy range. If in the genuinely 
Euclidean calculation one can reliably judge the accuracy achieved by considering 
the retained correction terms, an estimate of the accuracy for Minkowskean averages 
(integrals over a finite energy range) based on individual condensate terms is grossly 
misleading. To see that this is indeed the case it is sufficient to consider a model 
spectral density suggested in the last section of Ref. [TPl. This spectral density 



corresponds to an infinite series of equidistant poles, with one and the same residue, 
and it may be relevant in the large N c limit in the channel with one heavy quark 



2 



and a massless antiquark, 



1 00 1 
n mode i(£) = —ip(—E) + _ = -£ + Const , 

^ n=i ^ - n 

where ip is the logarithmic derivative of the T function. In order to mimic the 
standard QCD routine in the treatment of the spectral density we proceed as follows. 
One considers H{E) at Euclidean (negative) values of E, expands in 1/E and then 
takes the imaginary part of the expansion, 

ntsH-M-^ + g^^L (1, 

Im 11(E) = * - n -S(E) + * £ { ~ 1 , ) "'! l B " S l2 °' 1) (E) (2) 

Z n>0 \ Zn )- 

where B n stand for the Bernoulli numbers. The constant term in ImU(E) is an 
analog of the "perturbative" term; the rest is due to "non-perturbative" power 
corrections. Being extremely simple, the model contradicts no general requirements. 

It is not difficult to check that the Borel transform of this function, 511(e), 
considered at Euclidean (negative) values of the Borel parameter e, possesses the 
property we expect from the OPE-based analysis; namely, the exact result differs 
from the truncated series by a quantity of order of the last power term kept. At the 
same time the integrals 

Eo 

E n dE ImU 







calculated from the expansion (0) differ from the exact values by a large amount, 
~ 1/Eq, in spite of the fact that the power series for each given moment consists 
here of a finite number of terms. Thus, in the first moment (n = 0), if the first power 
term is retained, the second one and all others are zero, and one would expect the 
absolute accuracy following the line of reasoning of Ref . Jl3| . 



Of course, the model considered is relevant only in the limit N c = 00. For N c = 3 
a natural broadening of the resonances smears the spectral density and improves the 
accuracy of the Minkowskean calculations done with the truncated series. Still, the 
strength of the non-perturbative effects in the spectral density is represented not only 
by individual power terms of low dimensions but, also, by the asymptotic behavior 
of the power terms of high orders. The corresponding contribution is exponential in 
energy, 

AImTloce- CE ~ exp{-Ce c ' /as{E) } (3) 

and, therefore, is not seen in the truncated OPE expansions. Equation (|3p gives an 
estimate of a deviation from duality which shows up when one descends from the 
asymptotically high to lower energies. C and C in are constants, see Ref. ffM and 



the forthcoming publication [[HJ where the issue will be discussed in more detail. 
At the moment no reliable purely theoretical method exists that would allow one 
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to find the constant C. In other words, theoretical estimates of non-perturbative 
contributions in the Minkowskean quantities of the type of the total hadronic width 
of t or the e + e~ annihilation cross section at a given energy are rather vague, 
and they definitely do not have such a great accuracy as is required today in the 
problem of a s . One has to invoke phenomeno logical information, and this reverses 
the problem. If we accept that Aqcd lies in the vicinity of 200 MeV we have to 
conclude that about 20% of the pre-asymptotic term in the hadronic r width comes 
from non-perturbative effects (violations of duality) so that actually a s (M T ) « 0.27. 
At the moment this conclusion must be considered as perfectly legitimate. (By the 
preasymptotic term I mean R T — 3 where R T is defined in [jT3| .) Similar and even 
larger violation of duality was shown [16[] to take place in the inclusive semileptonic 
decays of the D mesons which are very close in mass to r's. 

What is usually done in the conference talks and review papers to lull the public 
opinion is averaging of two groups of data - low-energy and high-energy values of 
a s (Mz). Then the world average usually quoted is 0.117 ± 0.005; it lies only ~ 1.5 
standard deviations from either of them, and the contradiction is hidden under the 
rug. 

If the value of a s (Mz) ~ 0.11, as it stems from the low-energy data, what is the 
way out? The most placid solution of the problem would be reversing the trend of 
the Z peak experiments. Only three years ago the corresponding global fits used to 
yield numbers for the strong coupling constant which did not contradict the above 
value. Since then the result was steadily increasing. 

To bring the value of the strong coupling constant in line with the low-energy 
considerations one has to diminish the experimental number for the hadronic width 
of the Z by ~ 7 MeV. Surprisingly, this 7 MeV is the excess of the hidden beauty 
produced at Z, compared to the SM expectations, detected recently jl], As was 
noted in Refs. [[17], [l^] , if one allows for new physics in the Zbb vertex to take care 
of this 7 MeV excess in T(Z — > bb) one then solves the a s puzzle too. 

If no systematic bias is found (and experts say that this scenario is very unlikely) 
we are forced to look for physical explanations of the discrepancy. The SM global 
fits can be altered if there is a contribution due to new physics. As a matter of 
fact one of the explanations has been already proposed - light gluinos JITJ. Gluinos 
with masses of order of a few GeV change the rate of running of the strong coupling 
(it becomes slower), so that both numbers, the "low-energy" a s and the "high- 
energy" one become compatible with each other and compatible with the estimate 
Aqcd ~ 200 MeV. Many theorists, however, are reluctant to accept this scenario 
because of certain specific problems associated with the light gluinos. The niche for 
their existence is nearly closed experimentally. 

It seems more appealing to assume that new heavy particles (with mass > 100 
GeV ) generate, through loops, a correction to the Z-quark-antiquark vertices, en- 
hancing the hadronic decays of the Z. In order to convert a s [Mz) = 0.125 into 
a s (Mz) = 0.11 it is sufficient to ensure the enhancement of the hadronic width by 
~ 0.4%. The first idea that comes to one's mind is the fourth generation [2TJ. Then 
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the Z boson has a well-defined axial coupling to the doublet of quarks belonging 
to the fourth generation; let us denote them by T and B. Through the T and B 
loops the Z boson proceeds into a pair of gluons which are then coupled to light 
quarks, u, d, s, c, b. (One of the gluon propagators is contracted into a point due to 
the Z-boson quantum numbers). Interference of this graph with the tree Z-quark- 
antiquark vertices produces a correction to the hadronic decays of Z proportional to 
In mr/niBi so that the sign of the effect can be adjusted at will. The elegance of this 
mechanism becomes obvious if one takes into account the fact that the two-gluon 
intermediate state is weak isosinglet, so that the interference in the uu + dd channel 
cancels. The same happens in the ss + cc channel. The only surviving correction 
is in the bb channel (t is too heavy to appear in the intermediate state and cancel 
it). Thus, the fourth generation can naturally enhance the Z decays into bb, the 
mode where the current experimental data are known to disagree with theoretical 
expectations at the 2cr level [[I], |2J. If one could adjust the ratio In m^/m^ in such a 
way as to totally erase the disagreement in the bb then this bb enhancement would 



be sufficient to simultaneously solve the a s problem. Experts say, however J21 , |22 



that a large ratio mr/mg is ruled out by consistency of the SM radiative correc- 
tions to the masses and polarization operators of the Z and W bosons. Stretching 
all numbers to their extremes I found that the fourth generation can be responsible 
for at most 1.5 MeV in YiZ — > bb), instead of the desired 7. 

Leaving the fourth generation aside we can turn to superpartners. The light 
gluino scenario has been already mentioned. The heavy (virtual) gluino effects 
have been also discussed in the literature. The correction is generated by the Z 
coupling to squarks which then exchange a gluino and convert into quarks. Both 
squarks and gluino are assumed to lie in the 100 GeV ballpark. This mechanism 



was studied in [23], and later, even in more detail, in [24]. The squarks/gluino 
effect in the bb channel is specifically addressed in the works |25| . As it follows from 
these calculations, the sign of the gluino contribution is correct (i.e. it produces an 
enhancement) and, moreover, the effect can reach the desired 0.4% in the Z hadronic 
width provided that the gluino/squark masses are on the light side of the allowed 
mass domain. Thus, superpartners in loops can, in principle, solve both difficulties 
simultaneously. 

It should be noted, though, that if the superpartners are responsible for bringing 
a s (Mz) down to 0.11 the possibility of a straightforward Grand Unification within 
the Minimal Supersymmetric Standard Model (MSSM) is ruled out. Indeed, the 
simplest version of Grand Unification, with the squark and gluino masses in the 100 
GeV ballpark, implies ]26|] that a s (Mz) = 0.125 or larger. This seems to be an 
exciting observation defying the standard boring great desert scenarios [27|. 

Of course, one can say that the a s [Mz) problem - 0.125 versus 0.11 controversy - 
is only a 3<r effect. Being translated in the language of Aqcd the difference becomes 
quite drastic. Moreover, having Aqcd in the ballpark of 200 MeV is crucial for 
consistency of a very large number of Q CD-based calculations in the low-energy 
domain known to produce successful predictions which, seemingly, will not survive 
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if Aq CD > 500 MeV. In any case, the question is definitely ripe enough for an 
intensive public debate. The data on a s (Mz) from the Z peak, as they exist now, at 
the very least must be taken as a clear hint that new physics is around the corner. 



Whether it is supersymmetry or something else |17]] is hardly possible to decide 
at the moment. The supersymmetric explanation is good since (i) it goes through 
without twisting arms, i.e. it can naturally enhance the total hadronic and bb widths 
of the Z by the desired amount without spoiling the rest of the well-fit picture with 
the electroweak radiative corrections (see e.g. [pBfl); (ii) supersymmetry is anyway 
a popular element of the present-day theory as the only available mechanism which 
might explain the lightness of the Higgs particles in a natural way (assuming, of 
course, that relatively light Higgs particles do exist). A minimal lesson one has 
to draw is important for QCD practitioners. At the moment it seems reasonable 
to abstain from using a s {Mz) = 0.125 (and the corresponding value of Aqcd) as 
the best measurement of these key QCD parameters. If so, the problem of the 



semileptonic branching ratio deficit in the B mesons [^9] , which nearly disappeared 



after the corresponding theoretical formulae were evaluated with the large a s , 
resurfaces again. It would be interesting to check whether penguins generated by 
superpartners can ensure sufficient enhancement of the non-leptonic modes of B. 
Work in this direction has already begun, with quite encouraging results pi . The 



squark penguins can give rise to the chromomagnetic operators of the type sa^G^b 
with the coefficients less suppressed compared to the standard model expectations 



32 1 . Then the b — > s + gluon transition can be responsible for, say, 20% of the 
hadronic width of the B meson, thus eliminating any difficulties with Bt s i(B) pTf . 
Simultaneously, the expectation value of the charm multiplicity in the B decays 
goes down, which is also welcome. Further analysis is needed to check the overall 
consistency. It is necessary to verify, for instance, that the b — > sj rate is not 
enhanced beyond what is acceptable. Similar squark penguins can play a role in the 
AI = 1/2 rule in the strange particle decays. 

Another question to be considered is as follows. If the violation of duality at 
r is at the level of 20% of the pre-asymptotic term, what is to be expected in the 
inclusive B decays? 

All questions discussed above have been repeatedly considered in the literature 
previously - different aspects in separate publications. I combine them together. 
The only element which I add is my deep conviction that Aqcd can not be larger 
than ~ 200 MeV. The measurements of a s at the Z pole must be interpreted as a 
direct indication on new physics. Since convictions are very hard to formalize this 
letter should be viewed as an open invitation to further discussions among experts. 

Illuminating discussions with K. Hagiwara, H. Ohnishi, L. Roszkowski, M. Vir- 
chaux and especially P. Langacker, V. Novikov and M. Vysotsky are gratefully 
acknowledged. I would like to thank B. Holdom, P. Langacker and M. Vysotsky for 



pointing out to me Refs. |T7|, 113, [22|, E5f. This work was supported in part by DOE 



under the grant number DE-FG02-94ER40823. 
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